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Abstract We analyse a stochastic SIS epidemic amongst a finite population parti- 
tioned into households. Since the population is finite, the epidemic will eventually go 
extinct, t.e. have no more infectives in the population. We study the efi'ects of pop- 
ulation size and within household transmission upon the time to extinction. This is 
done through two approximations. The first approximation is suitable for all levels 
of within household transmission and is based upon an Ornstein-Uhlenbeck process 
approximation for the diseases fiuctuations about an endemic level relying on a large 
population. The second approximation is suitable for high levels of within household 
transmission and approximates the number of infectious households by a simple homo- 
geneously mixing SIS model with the households replaced by individuals. The analysis, 
supported by a simulation study, shows that the mean time to extinction is minimized 
by moderate levels of within household transmission. 

Keywords SIS epidemics • households model • time to extinction • Ornstein-Uhlenbeck 
process 



1 Introduction 

Epidemic models are widely used for increasing the understanding of infectious dis- 
ease dynamics {e.g. Anderson and May [1]) and for determining suitable preventive 
measures to reduce or ideally stop disease spread [e.g. Wallinga et al. [19], Ferguson 
et al. and Cauchemez et al. [^). In the present work we study a class of models 
for diseases that are endemic. More precisely we look at SIS-epidemic models (Kryscio 
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and Lefevre [13]) meaning that individuals are either susceptible (S) or infectious (1), 
and susceptible individuals might get infected and infectious, and after a time in the 
infectious state individuals recover and return to the susceptible state. The community 
is considered closed being of size A^, but an alternative interpretation of the model is 
that an individual that stops being infectious becomes immune for the rest of its life 
(or dies) but is "replaced" by a new susceptible individual when the infectious period 
stops. 

SIS-epidemic models can potentially lead to the disease establishing itself within 
the population and becoming endemic for a long period of time. However, eventually 
the disease goes extinct. It is therefore important to study properties of the model that 
determine whether the disease may become endemic, and if so what the endemic level 
is and for how long (on average) it will take for the disease to go extinct. In terms 
of prevention the answers to these questions can give guidance in determining how 
effective a given preventive measure is, for example: what proportion is necessary to 
vaccinate in order to make the disease go extinct promptly? 

In the present paper we address the question of the time until an endemic disease 
goes extinct, where we assume the community size A'^ is large and that we start in the 
endemic level. More precisely we study the situation where individuals are grouped 
into households and where it is assumed that transmission is much higher between 
individuals in the same household as compared with individuals in different households. 
We derive an approximation for the expected time to extinction and compare it with the 
case that there are no households (homogeneous mixing) for which asymptotic results 
for the time to extinction are already available (Andersson and Djehiche [3]). This 
paper is hence concerned with analysing how the time to extinction is affected when 
the assumption of homogeneous mixing of identical individuals is relaxed by admitting 
the presence of households. Similar questions have been analysed when relaxing the 
assumption of homogeneous mixing but in other "directions": Hagenars et al. Tl look 
at the case with an SIR epidemic where the community consists of a few large sub- 
communities assuming higher contact rates within the sub-communities, Lindholm [l^ 
studies an epidemic model where he investigates how the time to extinction is affected 
by individuals being of different types having different susceptibility and infectivity to 
the disease. 

Our main conclusion is that, assuming a fixed endemic level of infectives i, the time 
to extinction is smaller if there is moderate or high transmission within households as 
compared to the case of homogeneous mixing in the limit as A'^ — > oo. For high values of 
i our Ornstein-Uhlenbeck approximation suggests a monotone decrease of the expected 
time to extinction as a function of the degree of within-household transmission rate 
A^, whereas for small i-values the expected time to extinction first increases and then 
decreases as a function of A^. In the latter case the longest time to extinction is hence 
for the case with small (but positive!) Xj^. The Ornstein-Uhlenbeck approximation is an 
asymptotic result as N ~> oo using a normal approximation for the probability of being 
close to extinction. We introduce a second approximation which is appropriate when 
the within-household transmission rate, Xj^ , tends to infinity. This approximation shows 
that for fixed population size A'', the mean time to extinction grows polynomially with 
X]^ as — > oo with the rate of the polynomial growth depending upon the maximum 
household size. Thus there are two competing approximations leading in practice to 
the mean time to extinction being smallest for moderate levels of within household 
infection. We focus the analytical analysis on the case with households of size 2 with a 
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simulation study showing that qualitatively similar results hold for households having 
other and possible different household sizes. 

In Section 2 we define the SlS-household epidemic model and present some known 
results for the model (Ball [5], Neal [IS]). In Section 3 we derive our two approximations 
for the expected time to extinction and give a full discussion of the merits of the two 
approximations. In Section 4 we give some numerical examples and plot the expected 
time to disease extinction as a function of the degree of transmission within households, 
and compare these approximations with simulation results. The paper ends with a 
short discussion about the limitations of the present analysis, and a discussion about 
interesting related models/questions. 

2 The SIS household epidemic model 

2.1 Model definition 

We now define the Markovian household SIS epidemic model. This a special case of 
the more general household SIS epidemic model defined and analysed by Ball [5]. 

Consider a fixed community consisting of n households, all being of the same size 
k for simplicity of exposition (later we will primarily focus on the case k — 2). The 
community size is hence A'' — kn, and we assume that n and hence A'^ is large. Each 
individual is at any time either infectious or susceptible. Susceptible individuals recover 
(and become susceptible again) independently at the constant rate 7. While infectious 
an individual makes two types of contacts: the individual has global "close contacts" 
at rate Ac each time the contacted person is selected uniformly at random from 
the whole community, and the individual has local "close contacts" at rate with 
individuals belonging to the same household, here too the individual to be contacted 
is selected uniformly at random, but now among the individuals belonging to the same 
household as the person in question. By "close contact" we mean a contact that results 
in infection if the other person is susceptible - if the other person is infectious the 
contact is assumed to have no effect. The epidemic is initiated by some given initial 
condition specifying how many households have j infective individuals, j = 0, 1, . . . , fc 
infective individuals. 

Let Yj (t) denote the number of households having j infectives (and hence k — j sus- 
ceptibles) at t. The initial condition is hence specified by the vector {Yq{0), . . . , yfe(O)), 
for any t: Yj (t) — n since the total number of households is n. Let I{t) = jYj (t) 
denote the number of infectives (^infectious individual) at t, and similarly S{t) — 
— j)Yj{t) is the number of susceptible individuals. 

The process (Yo(t), . . . ,Y]^{t)) is Markovian and there is at most one infection or 
recovery occurring at a given time point, implying the possible jumps are that one 
component decreases by 1 and at the same time the component directly to the left 
(recovery) or directly to the right (infection) increases by 1. From the model, each 
individual recovers at rate 7, so the rate at which some household having i infectives 
at present has a recovery (implying that it changes to state i — 1) equals 7 • iYi{t). An 
infection can occur in two different ways. Either an individual is infected globally, or 
else locally from within the household. The rate of being infected globally is the same 
for all susceptible individuals: XQl{t)/N , each infective has close contact with a given 
individual at rate \g/N and there are I(t) infectives (in principle the division should 
be with A'^ — 1 since it is not possible to have contact with yourself, but when A'' is 
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large this is negligible) . Additionally, a susceptible in a household having i infectives at 
present gets infected locally at rate iX^Kk — 1) since each infective has close contact 
with a specific household member at rate Ai/(fc — 1). The overall rate at which some 
household having i infectives moves to the state i + 1 is hence 

(\Gl{t)/N + iXL/{k-l)){k-i). 
2.2 Properties of the SIS household epidemic model 

From now on we restrict ourselves to the case where all households are of size k — 2. 
The reason for this is to simplify notation and keeping the dimensions of vectors and 
matrices low - there is nothing harder in principle for larger household sizes. The SIS 
household epidemic model have been analysed previously be Ball [5] , Ghoshal et al. [10] 
and Neal ^5\, in particular studying how the model behaves when n is large. To this 
end we define Yi(t) = Y.i{t)/n, the proportion of households having i infectives and let 
I{t) = I(t)/N, the proportion of the population infectious. In Ball it is shown that, 
as n — > oo, the process {YQ{t), Yi {t),Y2{t)) converges in probability (uniformly on finite 
intervals) to the deterministic vector process {yo{t),yi{t),y2(.i)) defined as the solution 
to the differential equations 

y'f){t) = yyiit) - \Gi{t)2yo(t) 

y'l{i) = \Gi{i)'2yo{t) + 272/2 W - 7^1 (i) - (Ag^W + Wyi{t) 

y2 W = -27y2(t) + (AGi(t) + \L)yi{t), 

where i(t) — ^^ii/j(t)/2 (the fraction of infectives). The result relies on the initial 
conditions agreeing and having a positive fraction of infectives: (yo(0), yi(0), j/2(0)) ~ 
(yo, 2/1, J/2) as well as C?'o(0), ?l(0), ^2(0)) ^ (2/0,^1,2/2) as n ^ 00, such that J^j Vj = 
1 (the proportions sum up to 1) and '^jjyj/^ = j(0) > (there is a positive initial 
fraction of infectives). 

Ball [5] also derives a threshold parameter i?* = {XQ/'y)(l + \lH) such that 
the deterministic epidemic will go extinct (i.e. j/o(^) — > 1 as t 00) if i?* < 1, 
whereas it will reach an endemic level if > 1. In the latter case the endemic level 
(^0,^1,2/2) is given by Ball [S], Theorem 3.1. This gives that the endemic level of 
infectives i — ^jjyj/'i equals 

■ = 2A^7 7 + 

Ag(\/(Al + Ag)2 + 4Ai7 - (Ai + Ac)) Ag 

3 The time to extinction for the SIS household epidemic 

In the present section we derive a rough approximation for the time to extinction of the 
SIS household epidemic model defined in Section [2.11 applicable for all values of A^ and 
relying on n being large, and a more precise approximation suitable for large values 
of A^. Before doing this we present the results of Andersson and Djehiche [4] for the 
time to extinction of the homogeneous SIS epidemic model, i.e. the model above but 
without having household structure corresponding to the case where A^^ = 0. Without 
loss of generality we shall fix 7 = 1 which simply means that we measure time with 
the average infectious period as the base time unit. 
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3.1 The time to extinction for the homogeneous SIS model 

When there are no households ("the homogeneous case") the process is described by 
{S{t), I{t)), the number of susceptible and the number of infectives, but since S{t) + 
I{t) = N is fixed it is enough to keep track of I{t) say. When A'' is large and assuming 
R* — Xq > 1, I{t) will fiuctuate around the endemic level Ni — N{1 — 1/\q) but 
eventually make a big "excursion" down to I{t) — when the epidemic goes extinct. 
Andersson and Djehiche [4] prove that, as A'' tends to infinity, the time t^^^ until 
this happens is exponentially distributed with mean 

FC-r^^^h ~ , jV(lQgAG-l + l/AG) z^-, 

^^^^ ^ V^(Ag-i)2 ■ ^' 

The idea behind the proof is that the process will make many excursions from the 
endemic level eventually returning to the endemic level. Finally it will make a big 
excursion down to the absorbing state I(t) =0 when the epidemic goes extinct. It will 
hence make a geometric number of returns to the endemic level before going extinct, 
and in the proof they derive the expected duration of an excursion which returns to 
the endemic level and the (small) probability of making a big excursion to extinction. 



3.2 Approximation of the SIS household epidemic 

As mentioned in Section 12.21 it was shown by Ball [5] that when A'^ is large the SIS 
household epidemic may be approximated by a deterministic process. Using theory for 
Markov population processes (Ethier and Kurtz [8]) it is also possible to show a weak 
convergence result. Our process has three components, but because the total number 
of households X^^_g i^j(i) ~ n{= N/2) is fixed it is really two dimensional. This means 
it is enough to keep track of two components, e.g. Yq and Yi . Since our prime interest 
is the total number of infectives, it is natural to follow House and Keeling [12j and 
parameterise the model in terms of I{t) and Yi{t). Suppose the process is started from 
the endemic level, i.e. (/(O), Yi(0)) = {Ni, Nyi/2), and define the standardized process 

UrwJ-" Uiw/2-yi/2 

Then it follows from Ethier and Kurtz [8] that the vector process Z^^\t) with compo- 
nents m{t) andy/~)(t) converges weakly to Z, where Z is an Ornstein-Uhlenbeck 
process. The properties of the Ornstein-Uhlenbeck process stem from the system of 
differential equations ([T}, where we now have assumed that 7 = 1. These differential 
equations can in vector form be written as z' {€) — F{z{t)), where z{t) = {i{t),yi{t)/2). 
Using this notation, Z is defined by the drift matrix 

R = aF=f AG(l-2j(t))-l \l \ 

UG(l-2i(t)) + l-AGyi(t)/2-l-Ai,y ' 

and local covariance matrix 
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where f3i{t) = XGi{t){l - i{t) - yi{t)/2), P2{t) = {\Gi{t)yi{t)/2 + A^yi W), P:i{t) = 
2{i{t) — yi(t)/2) and P^^it) — yi{t) are the infinitesimal transition rates of an infection 
within an household with no infectives, an infection within an household with one 
infective, a recovery within a household with two infectives and a recovery within a 
household with one infective, respectively. Note that, if (j(0),yi(0)) = {i,yi), then for 
aU t, {i{t),yi{t)) = ilyi). 

This Ornstein-Uhlenbeck has a Gaussian stationary distribution with mean-zero 
and covariance matrix S defined by 

BE + EB^^~S. (6) 

This means that, for large t, our original process Z{t) will, conditional upon not having 
gone extinct, be approximately normal with mean vector {Ni, Nyi/2) and covariance 
matrix NS. 

Solving S is straightforward. However, the expression for S in terms of Xq and 
is not insightful as can be seen by the expression for En, the variance of the total 
proportion of infectives, given by, 

^ ('?2 - ft) J |r„^ 

" ((AL + AG)2 + 4Ai)<;3-?4 ' ^ ' 

where 

Ci = 6X1 + 22XgXI + 36Ai + 52X1 + 30XgXI + 76XgXI + 18X^x1 + 52Xlx% + 8X1 

+42XIXg + WX%Xl + 'iX^XL + 14A^Al + 2X% 
?2 = 2v'(AL + AG)2 + 4Ai(3Ai + SAgAI + 12Ai + 7X1x1 + IGXIXg + 8X1 + 2A^Ai 

+3ALAG + 6A^Ai + A^)) 
?3 = (5Ai + 12AGAi + 23Ai + 9AiA^ + 23Ai + 31AiAG + 13A^Ai + 2Xl + 2A^Ai 

+14AlAg + 2A?; + 2Ag) 
?4 = V(-^L + AG)2 + 4Ai(15A^Ai + 21A^Ai + 48AiA^ + 59Ai + 33Ai 

+5Ai + 22Ai + 2A^ + 2Ag + 68AGAi + ITAgA^ + GA^Ag 

+20A^Ai + SQAIAg + llA^Ai + 2A^Ai)). 

Similar expressions exist for the other components of S. For k > 2, the corresponding 
expression to ((Tjl is even more unwieldily. However, studying the behaviour of ((Tj in 
the limits as A/^ | and A^ — > oo is informative. 

Let cr?(i, Ai) = X'l^i, the variance for the total proportion of infectives in equilib- 
rium, explicitly stating the dependence upon i and A^^. (Note that Ag can be expressed 
as a function of i and A^.) Then for fixed i, crf{i, 0) = and Iim_!^^^oo '^fih Xl) = 
2i~^ . That is, for large A^ the variance of the proportion of infectives is approximately 
twice the corresponding variance in the homogeneous case. There is a simple explana- 
tion for this. In the homogeneous case, A^ = 0, we have N — 2n individuals. On the 
other hand, when A^ — > oo, the two members of a household are effectively paired to- 
gether, either both susceptible or both infectious. Thus the population in effect consists 
of n paired individuals resulting in the variance doubling. 

This leads onto the question of, whether or not, for fixed i, crf{i, Xj^) is monotoni- 
cally increasing in A^. Plots of rjf{i,X]^), fixing i and varying A/^, suggest that this is 
the case if i > 0.5, whilst, af{i,Xi^), is initially decreasing if i < 0.5. This is partially 
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confirmed by studying -;^^<^i{h^L)\^ _p wliicli is positive, equal to and negative 
when i > 0.5, i = 0.5 and i < 0.5, respectively. Note that for general choices of k, 
af{i,0) = i'^ and limx^->oo o-'f >^l) = ki~^ . 

If instead we consider fixed R* — Xci^ + ^ll is increasing as Xj^ increase, 
for all > 1. Also i is maximised at the extremes = and A^^ ^ oo, where 
1 = 1- 

The above analysis gives a good description of the endemic level, we now turn to 
the question of time to extinction and look to see how the above Ornstein-Uhlenbeck 
approximation can be used to assist in estimating this quantity. 



3.3 The time to extinction for the SIS household epidemic 

It seems hard to derive an explicit result corresponding to that of Andersson and Dje- 
hiche [3j for the household epidemic, the reason being that the process need not return 
to the endemic level {Ni, Nyi/2), and also because the trajectory down to extinction is 
not unique. It seems possible to derive a large deviation result but not to obtain a useful 
explicit expression for the time to extinction. Instead we have taken the approach first 
used by Nasell [T7] who applies it to the homogeneous SIR epidemic with demography, 
also used in Andersson and Britton We now present this approximation. 

Let Q denote the quasi-stationary distribution of the SIS household epidemic, i.e. 
Q is the stationary distribution of the process Y{t) = iYo{t),Yi{t),Y2{t)) conditioned 
on not having gone extinct. Hence, 

(y(t)|Y(0)~Q,/(t)>0)~Q. 

Starting in the quasi-stationary distribution it follows, because of the memoryless prop- 
erty, that the time to extinction Tq = inf{t; J(t) = 0|y(0) ~ Q} is exponentially dis- 
tributed with intensity equal to the probability of being one step away from extinction 
multiplied by the rate of moving into absorption. If we denote the quasi-stationary 
distribution Q = {qi^y^} = Pqilit) = i, Yi{t) = yi) we hence have that 

Tq ~ Exp(7gi,i) = Exp(gi,i), (8) 

since we have taken 7 = 1. Note that this is an exact result for any A'^. However, it 
remains to derive Since when I{t) = 1 we have by necessity that Yxify = 1 (and 
^o(^) = ^^"1)1 it is enough to look at the marginal distribution {qi} of the number 
of infectives. An approximation for qi is given by the normal approximation of the 
stationary Gaussian distribution derived in Section 13.21 This approximation consists 
of computing the normal density for I{t) (which is approximately normal with mean 
Ni and variance Naf {i, Xj^)) at the point 1 and conditioning on that I{t) > 0. This 
approximation is of course better in central parts of the distribution (around I{t) ~ Ni) 
but in the absence of a better approximation we use it also in the tail I{t) = 1. 

Using the expression for the endemic level (i, yi) and the variance (jf(i, Xj^) we get 
the following approximation for qi 

gi^P(/(t) = l|/(t)>0,Uarge) 

1 / 1 (l~Nif 

~ ■ exp — X 



V2^a,{i,XL) V 2 Naf{i,XL) 
1 / nP 



V2^a,{i,XL) V 2a2(i,AL) 



e^P . . . (9) 
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giving our first approximation for tlie time to extinction of the SIS household epidemic 
E{Tq) = - ^ V27NaS, Xl) exp ( \ . (10) 

The dominant term in our first approximation (IIOII is exp ( -^^ x -^r-L | . Thus 

for large A*' and a fixed value of i, the time to extinction is determined by af{i,X]^) 
with the smaller the value of crf{i, A^), the longer the time to extinction is expected 
to be. From the study of crf{i, A^) our first approximation hence suggests that E(Tq) 
is monotonically decreasing in A^ if i is fixed and larger than 0.5. If on the other hand 
we fix i < 0.5, E{Tq) first increases and then decreases with A^, and E{Tq) hence 
has a (local) maximum for a small but positive X^. It was noted in Doering et al. [7], 
that for the homogeneously mixing case the Ornstein-Uhlenbeck approximation is only 
reasonable for i?* = 1 + C/N^ , for C > and this is also likely to be the case for the 
household model. However, (|10p still proves to be useful in gaining an understanding 
of how Ai affects the extinction time. 

We now consider the case where A^ is large and derive a second more precise ap- 
proximation of E{Tq). As noted in Section [3.21 for large values of A^, it is unlikely 
that there will be just one infective in a household. Consider the infection of a sus- 
ceptible household making one member of the household infectious and assume for the 
moment that no further global infections are made with the household. Then after a 
waiting time of Exp(l -I- A^) either the second member of the household is infected 
(with probability X^ / (1 + Xi^)) or the infective recovers (with probability 1/(1 + A/,)). 
If the second member of the household becomes infected the time until one of these 
individual recovers is Exp(2). This recovery will shortly be followed by an infection 
(with probability Xi^/{1 + Xi^)) or a second recovery (with probability 1/(1 + A/^)). Let 
Ai,A2,... be independent and identically distributed according to A ~ Exp(2) and let 
Bq, Bi, . . . be independent and identically distributed according to B ~ Exp(l 4- A^). 
Let G ~ G'eom( 1/(1 + A^)) with support on Z"*". Let R be the time from infection of 
the household until it recovers and 5* be the sum of the total time infectious of the two 
individuals in the household from global infection until recovery. Then 



G 

R = Bo + ^iAj+Bj,, 

j=i 



where the sum is if G = and 

G G 



j=0 3 = 1 



Now X;'Lo ^3 ~ Exp(l) and 



^ r with probability 

_^ ^ ^ \ Exp(2/(1 + A^)) with probability 



Consequently, for large X]^, _R ~ _R ~ Exp(2/(1 -f- A^)) and whilst both members of 
the household are infectious the household is generating global infections at the points 
of a homogeneous Poisson point process with rate 2Xq. Since Xj^ is large, it is highly 
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unlikely that a global infectious contact (A^ is necessarily small) with a household 
with at least one infective in it will be with a susceptible individual. Therefore for 
large A^, the total number of infectious households (approximately the total number 
of infectives divided by 2) can be approximated by an SIS epidemic in a homogeneously 
mixing population of size N/2 with infection rate 2\q and recovery rate 2/(1 + 
From ([2]), this gives our second approximation for the time to extinction 



E{Tn)~ ^ + ^^ X W— (logi^.-l+l/i^.) 



where R* — Xg{^ + ^l)- Thus the time to extinction depends upon A^ and A^. This 
suggests that the mean time to extinction will approximately grow linearly in A/,, for 
fixed i, as A^ ^ oo. (For large Xj^, ~ 1/(1 — i).) 

This result can be extended to households where > 2. In general, for large 
values of A^^, the number of infectious households can be approximated by a ho- 
mogenously mixing SIS epidemic model with infection rate kXQ and recovery rate 
k } -, — , — T- For fixed A^, the mean extinction time will behave like CmX^~^ 

as Ai ^ CX3 for some Cjv > 0. Also lim;^^^^^^ cr^(j, A^) = ki~^. Thus the effect of 
A^ (households) on mean time to extinction is more marked for larger values of k. A 
similar result holds for unequal sized households. 

There is an apparent contradiction between the two approximations with the Ornstein- 
Uhlenbeck approximation (|10p stating that E{Tq) decreases in Xj^ for large enough A^ 
and the second approximation (|lip stating that E{Tq) increases in A^ as A^ oo. 
The explanation for this is that there are two asymptotic regimes considered: N ^ oo 
(Ornstein-Uhlenbeck approximation) and Xj^ ~* oo (second approximation). For fixed 
population size, A'^, there is a cross-over from the Ornstein-Uhlenbeck approximation to 
the second approximation as A^ increases. Hence this suggests that as X^ changes from 
to oo and when i < 0.5 (the most common situation) the mean time to extinction 
should first increase, then decrease and, when moving over to the second approxima- 
tion, eventually start increasing again. If we instead consider the situation where A*' 
increases, the transition between the two approximations occurs at increasing values of 
A^. Thus as A'^ ^ oo, the Ornstein-Uhlenbeck approximation dominates for all values 
of Xl. 

Equations ((2)1 and (llip give approximate mean times to extinction in the cases 
where A^ is close to and A^ is large, respectively. The question remains of estimating 
the time to extinction for moderate values of A^ since as noted by Doering et al. [7] 
the Ornstein-Uhlenbeck estimation of qi can be several orders of magnitude too small 
even for moderate A'^. However, the Ornstein-Uhlenbeck approximation can be used for 
a qualitative assessment of how the mean time to extinction changes when departing 
from the homogeneous case. That is, we can compare the Ornstein-Uhlenbeck approxi- 
mation for qi under the assumption of homogeneously mixing and for moderate values 
of X[^ for a given value of i. For example, for i < 0.5, does the mean extinction time in- 
crease for small values of A^ before decreasing as suggested by the Ornstein-Uhlenbeck 
approximation? Furthermore, what are the competing influences of A'^ and A ^ through 
the two approximations on the mean time to extinction? These questions are addressed 
in the following section. 
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Figure 1: Mean extinction times from 10000 simulations: Left figure N = 
N = 200. 



: 50 and right figure 



4 Numerical examples and simulations 

For the numerical examples and simulations we focus attention upon i = 0.2. In the ho- 
mogeneous case this corresponds to Aq = 1.25. Thus (7^(0.2, 0) = 5 and hm;^^^tx) o"^(0.2, A^;,) 
10. A simulation study involving iV = 50 and N = 200 and A^ = 0, 0.1, 0.2, . . . , 1, 2, . . . , 10 
was run, with for each set of parameter values the mean time to extinction estimated 
from 10000 simulations starting at the endemic level (see Figure 1). For — 50, the 
mean time to extinction was found to grow linearly in A^. The correlation between 
A^ and mean time to extinction was found to be 0.9992. The case Xj^ — 50 was also 
tested and found to satisfy the trend found for smaller values of A^. We conclude that 
the Ornstein-Uhlenbeck process (relying on N to be large) is not applicable for this 
case, and the second approximation is better for the whole range of A^. However, as 
N increased a different story emerged. First, for very small A^ E{Tq) seems to in- 
crease slightly and then decreases up until Xj^ ~ 2.5. Hence this part agrees with the 
behaviour suggested by the first approximation. After this, i.e. for Xj^ > 2.5, E{Tq) 
starts growing close to linearly with A/, as suggested by the second approximation pip . 



Similar results can be obtained where all households are of size > 2. In such cases, 
as A^ ^ oo, the mean extinction time for fixed A'' and i increases polynomially in A^. 
For moderate values of A^ and i < 0.5, the mean time to extinction is maximized by 
small (but positive) values of A^ as with households of size 2. This is illustrated in 
Figure 2 below where households of size fc = 5 are considered, and where the mean 
extinction time has been simulated for A = 50 and A = 200. The results are similar 
with those for k = 2, only more markedly. When A = 50 the second approximation is 
better for all A^ (suggesting polynomial growth of E{Tq) as a function of A/,) whereas 
when A = 200 the first approximation (Ornstein-Uhlenbeck), suggesting that E{Tq) 
should first increase and then decrease, works for small and moderate A^^, and after 
that, the second approximation suggesting polynomial growth in A^ performs better. 

In Figure 3, the ratio of the mean time to extinction with A = 200 divided by 
the mean time to extinction with A = 50 is plotted. By (|10p . this ratio should be 



N=50 



N=200 
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Figure 2: Mean extinction times from 10000 simulations in the case k = 5: Left figure N = 50 
and right figure N = 200. 



Ratio of mean extinction time N=200 v N=50 



Ratio of mean extinction time N=200 v N=50 



Figure 3: Ratio of mean extinction times from 10000 simulations for N ■ 
Left figure k = 2 and right figure k = 5 



200 and N = 50. 



approximately, 



V2Tr X 200(T(i, AL)exp(200r/2o-^(i, A^)) 
V27r X 50a{i,\L) exp(50iV2cr2(i, Al)) 



2exp(3/<7^(i,AL)). 



(12) 



Thus (|12p states we should expect to see the ratio between the extinction times ini- 
tially increase before decreasing as A^ goes from to infinity, corresponding to (T^(i, Xj^) 
initially decreasing before increasing as A^ ^ cxa. Figure 3 is consistent, at least quali- 
tatively, with the Ornstein-Uhlenbeck approximation. Similar results are observed with 
A'' = 500 with the Ornstein-Uhlenbeck approximation 'valid' for larger values of A^. 

For i > 0.5, we would expect to see that the extinction time decreases as A^ 
starts increasing from 0, in the limit as N ^ oo. Analysis of this case is not presented 
for two reasons. Firstly, in realistic situations i is likely to be small and at the very 
least less than 0.5. Secondly, the mean extinction time is very large even for small A'^. 
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For example, for i = 0.8 and N — 50, in the homogeneous mixing case E[r'^ 
4.18 X 10^'^ compared with 22.55 when i = 0.2. 



5 Discussion 

We present two approximations for obtaining the mean time to extinction in an SIS 
household epidemic model. Approximating for large values of A^^ the SIS household epi- 
demic model by a homogeneously mixing epidemic model with the households treated 
as individuals is informative and explains the trend observed for the mean time to 
extinction in the simulation study as X]^ increases. However, it is only really useful 
when is several orders of magnitude larger than Xq. The most interesting case is 
when Xq and are of the same order of magnitude. In that case we have resorted 
to the crude Ornstein-Uhlenbeck approximation to gain an understanding in the mean 
time to extinction. Although this approximation can severely over-estimate the mean 
time to extinction it does perform very well in giving a qualitative assessment of how 
the mean time to extinction varies with A^ . 

As mentioned in the introduction, better approximations are in principle available 
using large deviations, see for example, Shwartz and Weiss [18] . However, in practice 
it is difficult, if indeed possible, to get an explicit expression for the large deviations 
calculations. 

The SIS epidemic model is the simplest epidemic model which exhibits endemic be- 
haviour. It would be interesting to extend the above analysis to more realistic epidemic 
models with household structure. A prime example would be an SIR epidemic model 
with demography (births of susceptible individuals) extending the work of Andersson 
and Britton [5] to include household structure. 

It would also be interesting to consider the mean time to extinction in other SIS 
epidemic models with heterogeneous mixing of individuals. Examples include the great 
circle model (Neal [16]) and epidemics upon random graphs {e.g. Andersson [2]). 
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